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Abstract 

An elementary recursive model accounting for the quark spin in the fragmen- 
tation of a quark into mesons is presented. The quark spin degree of freedom is 
represented by a two-components spinor. Spin one meson can be included. The 
model produces Collins effect and jet handedness. The influence of the initial quark 
polarisation decays exponentially with the rank of the meson, at different rates for 
longitudinal and transverse polarisations. 



1 Introduction 

Present Monte-Carlo event generators of quark and gluon jets do not include the parton 
spin degree of freedom, therefore do not generate the Collins [IJ and jet handedness [2] 
effects. These are azimuthal asymmetries bearing on one, two or three hadrons, which can 
serve as quark polarimeters. However the asymmetries may strongly depend, in magnitude 
and sign, on the quark and hadron flavors and on the transverse momenta and scaled 
longitudinal momenta z of these hadrons. Therefore a good knowledge of this dependence 
is needed for parton polarimetry. Due to the large number of kinematical variables, a 
hadronisation model which takes spin into account is urgently needed as a guide. 

The semi-classical Lund ^Pq mechanism [3J , grafted on the string model, can generate 
a Collins effect [4J, but not jet-handedness. Here we propose a fully quantum model of 
spinning quark fragmentation, based on the multiperipheral model. It reproduces the 
results of the ^Pq mechanism and also contains the jet-handedness effect. 



2 Some recalls about quark fragmentation 

Figure 1 describes the creation of a quark "go" and an antiquark in e"'"e~ annihilation 
or decay, followed by the hadronisation, 

qo + q-i ^ hi + h2... + . (1) 

Looking from rigth to left, one sees it as the recursive process (see [5] and ref. 4 of [0]), 

go = go ^1 + Qi 4-momenta : ko = pi + ki , 

gi ^ h2 + q2 ki=p2 + k2 , 

qN-i hN + qN kN-i = Pn + kN ■ 

1 presented at XIII*'^ WORKSHOP ON HIGH ENERGY SPIN PHYSICS, DSPIN-2009, Dubna, 
Russia, Sept. 1-5, 2009 ; http://theor.jinr.ru/~spin/2009/ 
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Figure 1. Electroweak boson qq mesons. 

gjv = Q-i is a "quark propagating backward in time" and k]\r = —k{q^i). 
Kinematical notations : 

ko = fc(go) and are in the +z and — z directions respectively. For a quark, t„ = knx- 

For a 4- vector, = a° ± and = (a^, a^). We denote by a tilde the c?ua/ transverse 
vector = z X ax = (— a^, a^)- 

In Monte-Carlo simulations, the kn are generated according to the splitting distribution 

dW{ g„_i ^ hn + Qn) = /n(Cri, ^I^pIt^ ) ^^Cn C^^*n , Cn = Pt/K-l ■ 

In particular the symmetric Lund splitting function |3j, 

/„ oc -1 (1 - D exp [-6 {ml + p^^)/C„] , (3) 

inspired by the string model, fulfills the requirement of forward-backward equivalence. 

On can also consider [B] the upper part of Fig.l as a multiperipheral [7j diagram 
with the Feynman amplitude 

M-qo+q-x-^ln-.+hN = ^{k-l^ S^i) Tqj^^hj^^qj^_^{kN, kjy^i) Aq^_j (A;jv_i) ■ ■ ■ 

(/C2, ki) Aq^{ki) Tq^^huqoiki, ^o) ^(/co, ^o) ■ (4) 

So and S'_i are the polarisation vectors of the intial quark and antiquark. S'^ = 1, Sz = 
helicity, St = transversity. F and A are vertex functions and propagators which depend 
on the quark momenta and flavors. Note that Fig.l is a loop diagram : fco is an integration 
variable, therefore the "jet axis" is not really defined. Furthermore, in Zq or 7* decay, 
the spins go and g_i are entangled so that one cannot define Sq and S^i separately. 

Collins and jet-handedness effects. Let us first assume that the jet axis (quark 
direction) is well determined : 

- the Collins effect in g*— /i+X, is an asymmetry in sm.[ip[S)—(p{K)\ for a transversely 
polarized quark. The fragmentation function reads 

F{z,Pt]St) = Fo{z,pI) {1 + AtSt.Pt/\Pt\) (Pt = ^ x Pt) ■ (5) 
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Figure 2. String decaying into pseudoscalar mesons. 
At = At{z,Pt) G [—1, +1] is the Collins analysing power. 

- jet handedness [2\, in q ^ h + h' + X , is an asymmetry in sm[(p{h) — ^p{h')] proportional 
to the quark helicity. The 2-particle longitudinaly polarised fragmentation function is 

\ \Pt'Pt\/ 

Al = Al{z,p^, z' ,p'rp^ ,Px-Pt) G [—1, +1] is the handedness analysing power. Pi^ ■ p't is 
the same as z ■ {pr^ x p'r^)- 

If the jet axis is not well determined, an additional fast hadron, h' or h" is needed. 
The z axis is taken along P = {p+p') (Collins) or P = {p+p' +p") (handedness). In this 
way, we define the 2-particle relative Collins effect (also called interference fragmentation) 
and the three-particle jet handedness, which corresponds to the original definition of [2]. 

The Lund ^Pq mechanism . Figure 2 depicts the decay of the initial massive string 
accompagnied with the creation of a qq pairs. Forgetting transverse oscillations of the 
initial string, the transverse hadron momenta come from the internal orbital motions of 
the pairs. After a tunnel effect the q and g of a pair become on-shell and their relative 
position r = r{q) — r{q) is along — z. The pair is assumed to be in the ^Pq state, which 
has the vacuum quantum number. The relative momentum k = k{q) = —k{q) and the 
orbital angular momentum L = r x k are such that z ■ [k^ x L] < 0. In the ^Pq state 
(sq) = {sq) = —{L/2). As a result, the transverse spins of q and q are correlated to their 
transverse momenta : 

(z- [fcT(q) X Sq]) > 0, (z- [fcT(q) X Sq]) < 0. (7) 

The correlation can be transmitted to a baryon. Then ( z ■ [prj, x sb] ) has the sign of 
(Sq.s^). This can explain transverse spin asymmetries in hyperon production [5]. 

Application to the Collins effect [4j. In Fig. 2, go is polarised along the direction +y 
toward the reader and hi is a pseudoscalar meson, for which (■s(go)) = "(•^(gi))- Then 
gi and gi are polarised along — y and, according to ([7]), ^^(gi) = Prihi) is in the +x 
direction. This provides a model for the Collins effect. Fig. 2 also indicates that, for a 
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sequence of pseudoscalar mesons, the Colhns asymmetries are of ahernate sides. Besides, 
Qn-i and Qn go on the same side, which enhances the asymmetry. It may explain why vr" 
from M-quarks have a strong Collins analysing power. Note that this effect also enhances 
(p^), independently of the go polarisation. 

3 A simplified multiperipheral quark model 

In Eq.(jl]), let us replace Dirac spinors by Pauli spinors. A minimal model, restricted to 
the direct emission of pseudoscalar mesons, is built with the following prescriptions : 

1) replace u{kQ, Sq) and v{k_i, S^i) = —u{kq_-^, —Sq_J 75 by the Pauli spinors x('S'o) and 

2) assume no momentum dependence of F , 

3) replace 75 by , 

4) replace the usual pole (A;^ — m^)^^ of Aq{k) by the {ki, kr) separable form 

Dq{k) = Qgik+k-) exp(-5tV2) , (8) 

5) replace the usual numerator rriq + 'j.k by fiq{k~^k~ ,t^) + icr.t . 

These prescriptions respect the invariance under the following transformations : 

- rotation about the 2;-axis, 

- Lorentz transformations along the z-axis (longitudinal boost), 

- mirror reflection about any plane containing the z-axis (parity), 

- forward-backward equivalence. 

The jet axis being fixed, full Lorentz invariance is not required, whence the separate 
dependences of Dq and Hq in k~^k~ and t^. In item 5), fi+icr.t is reminiscent of the meson- 
baryon scattering amplitude f{s,t)+ig{s,t) cr.{pxp'). Single-spin effects are obtained for 
Sm /i 7^ 0. The choice of putting the spin dependence in the propagators rather than in 
the vertices is inspired by the ^Pq mechanism : in both models the polarisation germinates 
in the quark line between two hadrons. 

For a fast investigation of the model, we make the further approximations : 

- Neglect the influence of the antiquark flavor and polarisation in the quark fragmentation 
region. This is allowed at large invariant go + Q-i mass. 

- Discard the interference diagrams. For a given final state, the rank ordering of hadrons in 
the multiperipheral diagram is not unique and differently ordered diagrams can interfere. 
This interference (and the resulting Bose-Einstein correlations) will be neglected. 

- Disentangle k^ and k^. We will assume that ^q{k~^k~ ,1"^) is constant or a function of 
only. Thus we have no more "dynamical" correlation between longitudinal and transverse 
momenta. However there remains a "kinematical" correlation coming from the mass shell 
constraint 

{kn-l - = {k+_, - k+){k-^^ - k-) - [tn-l - = m\hn) . (9) 

In the following we will ignore the (t„-i — t„)^ term. This approximation is drastic for 
pion emission because (f^) > m^. We only use it here for a qualitative investigation of 
the spin effects allowed by the multiperipheral model. Thanks to it, the t's become fully 
decoupled from the k^ and kinematically decorrelated between themselves. They remain 
correlated only via the quark spin. 
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Py-distibutions in the quark fragmentation region. With the above approxima- 
tions we can treat the process ([2]), at least in pj,-space, hke a cascade decay of unstable 
particles, which has no constraint coming from the future. The joint distibution of 
the n first mesons is proportional to 

/(PiT,P2T, -Put) = expi-Btl -Btj ... - Btl) Tr !^M^2...n i±|^ Ml^.^^j , (10) 
with 

Mi2...n = M„ ■ ■ ■ M2 Ml , Mr = {fXr + io-.ir)a, . (11) 
3.1 Applications to azimuthal asymmetries 

In this section we will calculate azimuthal asymmetries for particles of definite ranks. 
For comparison with experiments, one should mix the contributions of different rank 
assignments. For simplicity we take a unique and constant fi for all quark flavors. 

First-rank Collins effect. Applying fllOtllip for n = 1 gives 

I{p,t) =exp{-Btl) + tl-2^m{fi)ii.S) , (12) 
with ti = —PiT- For complex /j, one has a Collins asymmetry (cf Eqj5]) with 

+Pl,T 

If ^m{fi) > it has the same sign as predicted by the ^Pq mechanism. 

Joint Pt spectrum of hi and Applying (11 QUI II) for n = 2 one obtains 

HPit,P2t) = exp{-Btl-Btl){i\fi\^ + tl){\fi\^ + tl)-Ati.t2Qm\ix) 

+ 2 55m(/i)5.ti(2ti.t2- l/xp-ta) 
+ 2 53m(/i)5.t2(|/i|'-*?) 

- 2S3m(/i2)5.(ti xt2)}, (14) 

with ti = —PiT, ^2 = ~{PiT + P2t) S ■ {ti X ±2) = 5*2 Pit ■ p^j,. 
The last line contains jet handedness (cf Eql6]), of analysing power 



— 2 3m(/i^) \piT X P2t\ 

(|/i|2 + t?)(|/iP + ti)-4ti.t2 W(/i) 



The second line contains the Collins asymmetry of hi. Both 2^^^^ and 3'"'^ lines contribute 
to the /i2 one after integration over ti, and to the relative 2-particle Collins asymmetry, 
which bears on 

ri2 = ■ = ■ t2-ti. (16) 

Zl + Z2 Zi + Z2 

Note that Collins and jet-handedness asymmetries are not maximum for the same value 
of arg(/i). This is related to the positivity [8] constraint 

^UPit, P2t) + Al{piT, P2t) < 1 ■ (17) 
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3.2 Evolution of the polarisation of the cascading quark 

Let us first assume that ti, t2, ■■■ tn are fixed and consider the spin density matrix 
Pn = {1 + Sn-cr)/2 of Qn at the {n + 1)**^ vertex : 

Pn = Rnl Tr{i?„} , = Mi2...„ Ml2...„ . (18) 

If po is a pure state (det po = 0), then p„ is also a pure state ; no information is lost. 

Let us now integrate over ti, t2, ... tn (equivalently over p^^, ...p,^j.). It leads to a 
loss of information. The spin density matrix of g„ becomes 

Pn = /?„/Tr{i?J, Rn = J d%...J dHn Mi2...n ^ ^ ^"'"^ MI^...^ • (19) 

Rn and Rn obey the recursion relations 

Rn = M„ Rn-, Ml, Rn = J dHn Mn Rn-i Ml . (20) 

At fixed t's, the left equation gives (setting p = p' + ip") : 

ir A'-I/^P -2\t\p'\ . 

5„= -<^2p"^„ + 7^[z,y^n] -\p?-t^ 7^[z,-V9n]Sn_lL (21) 

V 2|t|p' |pp-t7 ^ 

with C = Tr{_R„} = |pp + — 2p"tn ■ Sn-i. The rotation 7l[z,ipn\ about z brings x 
along tn. Iterating ( |T2l) . where we replace {S', ti} by {Sn-i, tn}, and ( 12T|) . we generate 
the successive transverse momenta with the Monte-Carlo method. From ( 12T|) we learn : 

- if Q'mp 7^ 0, the inhomogeneous term in p"tn is a source (or sink) of transverse 
polarisation : one can have SnT 7^ even with Sn~i = 0. 

- helicity is partly converted into transversity along tn and vice-versa. 

The last fact explains the mechanism of jet handedness in this model : first, the helicity 
5*^0 is partly converted into Sit parallel to Pirp, then Sit produces a Collins asymmetry 
for /i2 in the plane perpendicular to Pirp. 

Let us now consider the i-integrated density matrix. The right equation in ( 120|) gives 

Sn,z = DiiSn-l^z , Sn,T = Dtt Sn-l,T ', Dtt ^ [~1; +1] ; (22) 

(St) = / ^^P(-^*') C^-H / / ^^P(-^*') d^l' + • (23) 

Analytical values : Dll = (^ - l)/(^ + 1) and Dtt = -^/(^ + 1) with ^ = B\p\^. The 
geometrical decays of \Sn,z\ and \Sn,T\ along the quark chain occur at different speeds. 
They are similar to the decays of charge and strangeness correlations. Dn and Dtt 
saturate a Soffer-type [8] positivity condition 

2\Dtt\ <1 + Dll. (24) 

Indeed, 2Dtt = — 1 — -Dll- This is due to the zero spin of hn (compare with text after 
Eq.(4.87) of [8J). The negative value of Dtt leads to Collins asymmetries of alternate 
signs, in accordance with the ^Pq mechanism. It comes from the vertex for pseudoscalar 
mesons. For scalar mesons we replace by 1. In this case Dtt is positive, g„_i and g„ 
tend towards opposite sides and the Collins effect is small, except for hi. This is also the 
prediction of the ^Pq mechanism. 
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4 Inclusion of spin-1 mesons 



For a J = 1 vector meson and the associated self-conjugate multiplet, the "minimal" 
emission vertex written with Pauli matrices is 

T = GlV:1 + Gt (t.V*tCT,, (25) 

where V is the vector amplitude of the meson normalised to V.V* = 1. It is obtained 
from the relativistic 4-vector first by a longitudinal boost which brings the hadron at 
Pz = 0, then a transverse boost which brings the hadron at rest. 

For a J^*" = axial meson of amplitude A, the "minimal" emission vertex is 

T = Gt(t.A*t. (26) 

It differs by a matrix from the second term of f l25|l . A term of the form A* cr^ with 
constant is not allowed by the forward-backward equivalence. 

Let us treat the case where the l^*-rank particle is a meson and fix the momenta 
p{7i~^) and p(vr°) of the decay pions. Then oc p{tt~*^) —p{7T^), which is real, corresponding 
to a linear polarisation. Replacing the coupling of ( fTTj) by (!25|) we obtain 



I{Pt,V) = exp(-5t2) |G^|2x 

{ (|a| + Vl) (l/ip + t^) - AVr.t K 3m(a) 3m(/.) 
+ 2'^m{fi) \a\^V^ S.t 

+ 2'^m{fi){VT.iVT.S + VT.tVT-S) 

+ 4 3fJe(a) Qm{fx) 14 Vt-I } , (27) 



with t = ti = —prp[p^) and a = Gl/Gt- Let us comment this formula : 

• The 2"^^ line is for unpolarised quark. It gives some tensor polarisation. 

• The 3'''^ line is a Collins effect for the p"*" as a whole, opposite to the pion one (com- 
pare with 0121) and only for longitudinal linear polarisation, in accordance with the ^Pq 
mechanism. For (V^) = 1/3 (unpolarized p"*") and a = 1 one recovers the Czyzewski 
prediction [9] ^^(leading p)/^^ (leading vr) = —1/3. 

• The 4*^ line gives an oblique polarisation in the plane perpendicular to St corresponding 
to hi or hiLT of [inilll]. After p"*" decay, it becomes a relative 7r+ — vr" Collins effect. 

• The 5**^ line is a new type of asymmetry, in sin[2(y9(V) — ip{t) — ^p{S)]. 

• The last line also gives an oblique polarisation, but in the plane perpendicular to Pxip^)- 
After p"*" decay, it becomes jet-handedness. Indeed, ignoring an effect of transverse boost, 
we have Vt oc pj.(7r+) — pj,{iT^), therefore Vx-i oc —prp{'K~^) x Pj,(tc'^). 



5 Conclusion 

For the direct fragmentation of a transversely polarised quark into pseudo-scalar mesons, 
the model we have presented has essentially one free complex parameter p and reproduces 
the results of the semi-classical ^Pq mechanism : large asymmetry for the 2"'^-rank meson, 
Collins asymmetries of alternate sides for the subsequent mesons. In addition, it possesses 
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a jet-handedness asymmetry, generated in two steps : partial transformation of helicity 
into transversity, then Collins effect. 

We have also considered the inclusion of spin-1 mesons. When longitudinally polarised, 
a leading vector meson has a Collins asymmetry opposite to that of a pseudoscalar, as 
also expected from the ^Pq mechanism. The decay pions of a p meson exhibit a relative 
Collins effects as well as jet-handedness. These effects are associated to oblique hnear 
polarisations of the p meson. The fact that two pions coming from a p show the same 
spin effects as two successive " direct" pions is reminiscent of duality. 

Even for unpolarised initial quarks, the spin degree of freedom of the cascading quark 
has to be considered. It enhances the (p^) of the pseudoscalar mesons compared to scalar 
and longitudinal vector mesons. 

A next task for building a realistic Monte-Carlo generator with quark spin is to take 
into account the {tn-i — tnY term in ([9]). One must also be aware that there exist other 
mechanisms of spin asymmetries in jets. For example the Collins effect can be generated 
by the interference between direct emission and the emission via a resonance |12] . 



References 

J. Collins, Nucl. Phys. B 396 (1993) 161. 

A. V. Efremov, L. Mankiewicz, N.A. Tornqvist, Phy. Lett. B 284 (1992) 394. 

B. Andersson, G. Gustafson, G. Ingelman and T. Sjostrand, Phys. Rep. 97 (1983) 31. 
X. Artru, J. Czyzewski and H. Yabuki, Zeit. Phys. C 73 (1997) 527. 
A. Krzywicki and B. Petersson, Phys. Rev. D 6 (1972) 924; 
J. Finkelstein and R. D. Peccei, Phys. Rev. D 6 (1972) 2606. 
X. Artru, Phys. Rev. D 29 (1984) 840. 

D. Amati, A. Stanghellini and S. Fubini, Nuov. Cim. 26 (1962) 896. 
X. Artru et al, Phys. Rep. 470 (2009) 1. 
J. Czyzewski, Acta Phys. Polon. 27 (1996) 1759. 
Xiangdong Ji, Phys. Rev. 49 (1994) 114. 

A. Bacchetta and P. J. Mulders, Phys. Rev. D 62 (2001) 114004. 
J. C. Collins and G. A. Ladinsky (1994) arXiv : |hep-ph/94H444l 



8 



